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Abstract. 

An anisotropic t — J — U model Hamiltonian is used to model electron behaviour in 
quasi-2d materials in the dilute limit, and as a highly simplified representation of the 
weakly coupled Cu02 planes of the high-T^ cuprates we model the very poor out-of- 
plane conductivity via the complete suppression of interplanar hopping. However, we 
do include the very weak interplanar superexchange, and are thus considering a model 
of exchange-coupled planes. For an isotropic three-dimensional system in the dilute 
limit, we find that the formation of two-particle bound states requires Jc/t > 5.9. 
Also, it is known that Jc/t = 2 for a 2d square lattice. However, for our model of 
exchange-coupled planes any infinitesimal interplanar exchange (Jc — 0) is adequate 
to form bound states. 
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1. Motivation and Introduction 

The discovery of high-temperature superconductivity in copper-based transition metal 
oxides was made by Alex Miiller and Georg Bednorz in 1986 (see [1] for more 
information). They found that the copper oxide compound {La, Ba)2Cu04^ became 
a superconductor at Tc ~ 30K, which was substantially higher than for other any other 
compounds known at that time. Presently, one can find transition temperatures at 
ambient pressure close to 138K for other cuprate-based systems. The mechanism for 
this novel behaviour is still a subject of spirited debate, but one idea that has been put 
forward repeatedly is the Heisenberg superexchange between neighbouring Cu sites, a 
biproduct of the strong repulsive electronic (Hubbard-type) correlations that are present 
on the transition metal sites. 

The structure of all of the cuprate materials is similar, in that Cu02 planes are 
stacked one on top of another to produce a quasi-2(i crystallographic arrangement. 
This characterization of these materials is supported by their very poor interplanar 
conducting behaviour [2j, at least in the weakly doped regime [H Hj. The so-called 
c-axis puzzle has been studied extensively [5], and its relation to the superconductivity 
has been discussed [6]. Further, given the dependence of Tc on the number of Cu02 
planes per unit cell [7], the possibility of an interplanar pairing mechanism cannot be 
ignored, and indeed previous work has shown P[9], within various approximations, that 
any interplanar interaction increases T^. 

In this report we examine the two-electron problem {viz. the dilute limit) in the 
highly simplified situation of zero interplanar hopping. That is, our model is meant to be 
a very rough approximation to the extremely low out-of-plane conductivity, but leaving 
the residual interplanar superexchange found in the cuprates, the latter produced by 
strong electronic correlations. We refer to this situation as exchange- coupled planes. 
We solve for the conditions under which a two-particle bound state is formed, and in 
particular determine the minimum value of superexchange {Jc/t) for which bound states 
appear. Our results make evident the potential benefit of having electrons confined to 
individual planes that can only move via the spin-exchange {^S^S~) process. To be 
specific, we find that while for a two-dimensional plane {Jc/t = 2) or an isotropic three- 
dimensional system {Jc/t ~ 5.9) one requires a superexchange much larger than that 
found experimentally {J/t ~ 0.3), for exchange-coupled planes one requires only an 
infinitesimal {Jc = 0) interaction. Therefore, even the very small out-of-plane exchange 
coupling (estimated to be roughly 10~^ of the in-plane exchange) would suffice to form 
bound states. We note that similar results were found previously for exchange-coupled 
chains |10j, emphasizing the potential importance of such electronic confinement. 
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2. Formalism 

2.1. Model Hamiltonian 

We consider tlie so-called t — J — U or Heisenberg-Hubbard model Hamiltonian, given 
by 

In this Hamiltonian i,j label the sites of the lattice, the notation {i,j) denotes 
neighbouring lattices sites in any of the x, y ot z directions, each near neighbour pair is 
counted once only, and Ci^„/cl ^/rii^^ denotes the annihilation/creation/number operator 
for an electron at site i with spin a. We will consider the dilute limit, and in particular 
only two electrons, and therefore the frustrating geometry of some the high-T^ cuprates 
(such as occurs for a body-centred tetragonal structure) can be ignored. Therefore, we 
consider a simple tetragonal lattice structure, which reduces to a simple cubic structure 
for isotropic hopping and exchange - see below. 

For the energy parameters in the Hamiltonian we consider the following situations; 
a schematic of these parameters is shown in figure [H First, the hopping frequency is 
restricted to be near neighbour only, and in-plane and out-of-plane hopping frequencies 
are allowed to be different: 

tij = t|| when i,j are near neighbours in either the x oi y directions 
= t± when i,j are near neighbours in the z direction (2) 
Similarly for the superexchange integral: 

"^iJ ~ "^W when i,j are near neighbours in either the x or y directions 
= Jj_ when i,j are near neighbours in the z direction (3) 

The Hubbard on-site repulsion energy is parametrized by U. 

Strong coupling expansions of the Hubbard model show that to lowest order in t/U 
one has J = 4t'^/U, and in the t — J — U model this exchange interaction is included 
explicitly. Note, however, that this model Hamiltonian does not involve projection 
operators that project away doubly occupied sites, as one considers in t — J model 
studies [11]. Instead, we take the limit ?7 — *• cxo to represent the high cost of doubly 
occupied sites, similar to what has been done in previous studies [TOl [12] • That is, J 
is the effective interaction that is created by some finite value of U [13], and we then 
reduce the Hilbert space by taking the U ^ oo limit leaving J nonzero. 

As we have discussed in the introduction, in this paper we compare the bound 
states for a two-electron system, described by the Hamiltonian of ([1]), when interplanar 
hopping does not occur (that is tz = 0), to the situation in each of the hopping and 
superexchange energies in ([2]) and ([3]) are the same in all three directions. 

Consider the two-electron, 'S'^^^^j = state 

|^) = E0(m,n)aLXalO) (4) 

m,n 



Bound state formation in exchange- coupled planes 




± 




Figure 1. (Colour online) A schematic representation of the various terms that appear 
in the t — J — U model Hamiltonian. The electron hopping t, both in-plane and 
inter-plane {tj_), spin exchange J, again for both in-plane (J||) and inter-plane (J_l), 
and the on-site Coulomb repulsion U. 



created from the empty lattice, |0). With the added symmetry (f){m,n) = (j){n,m), {tp) 
will be a singlet state. We use the following for the Fourier coefficients for (f){i,j): 



Introducing ki = ^ -|- q and k2 = -y — q, we restrict our attention to zero centre-of-mass 



2 ' ^ ""^ 2 

momentum states, and therefore solve for solutions in terms of 



(6) 



0(ki,k2) = 0(q, -q) = 0(q) 
We make the following definitions: 

e(k) = - 2{t|| [cos(k ■ ic) + cos(k ■ y)] + t^ cos(k ■ z)} (7) 
J(k) = 2{ J|| [cos(k ■ x) + cos(k ■ y)] + cos(k • z)} (8) 
and find [121 that H\iIj) = E\tp) is equivalent to solving the following equation: 

, f Ep0(p)-^Ep^(q-p)0(p) 
'^(^^ = 

To simplify ([9]) it is convenient to make the following definitions (with Xi and xj 
being any of x, y or z): 



(10) 



M ^ 



N ^ E + 4(t|| [cos(q ■ x) + cos(q ■ y)] + t± cos(q ■ z)) 



:iii 



cos(q ■ Xi) 



iV ^ E 4(t|| [cos(q ■ x) cos(q ■ y)] ■\- 1± cos(q ■ z)) 
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AT ^ E + 4(t|| [cos(q ■ x) + cos(q ■ y)] + cos(q ■ z)) 

By substituting these definitions into iQ one straightforwardly obtains a system of 
equations that one can solve for the bound states of the problem: 

Co = UIoCo - 2J\\{hCx + lyCy) - 2J^/,a (12) 

Cx = UIxCq — 2 J|| (/j-a-Ca; + IxyCy) — 2J±IxzCz 

Cy = UlyCo 2J^\(^IxyCx + lyyCy) 2J±IyzCz 

Cz = UIzCq — 2J\\(IxzCx + lyzCy) — 2Jj_IzzCz 

3. Results 

We analyzed the above system of equations for two geometries of lattices: infinite in all 
directions (referred to as the 3d case), as well as for Lx = Ly = oo and Lz = 2 (referred 
to as the two plane problem). For the 3d case, where Lx, Ly, Lz oo, the sums in ffTOl) 
and (fTTl) become 



-E 



dqx / dqy / dqz (13) 



N ^ (27r)^ 

and similarly for an isotropic 2d lattice. Also, for Lx, Ly oo in the two plane problem 
{Lz = 2) the sums become 

For the t — J — U model with isotropic hopping {t) and exchange (J) in dimension 
d, the minimum (total) energy for two noninteracting electrons is -Adt, and therefore 
the energy of a bound state must be less than this energy. For convenience, we define a 
scaled energy A by 

A = ^ (15) 
At ^ ' 

One of our results will be related to the simpler and lower-dimensional isotropic 

single-plane problem, and thus here we outline a method of solving such a problem. 

Eliminating the z dependence and using the equivalence of the x and y directions, (IT^ 

reduces to 

Co = Uh,2dCo-^JIx,2dCx (16) 
Cx = UIx,2dCo — 2JIxx,2dCx — 2JIxy^2dCx 

One can further simply these equations using the identities 



h2d =1-^ (17) 

Ixx,2d = ~ Ixy,2d ~ AIx,2d 



Bound state formation in exchange- coupled planes 



6 



and one can integrate Io^2d to obtain 

1 



where K{x) is the complete elliptic integral of the first kind. Thus /o,2d is a symmetric 
function of A, which is defined for |y4| > 2 and diverges as — >■ 2+. Now we eliminate 
Co and from the equation set (one can show that no solutions exist for either of these 
being zero - see below for 3d), and then take the strong Coulombic repulsion limit of 
U oo, and then solve for J. One finds 

J=^l%^ (19) 

Noting the dependence of /o,2d on energy, this equation (and we will derive similar 
equations for 3d below) then solves for the bound state energy in terms of J. 
Alternatively, one can use this equation to solve for the minimum value of J/t that 
is required for a bound state to form. That is, by using ( |T8i) one finds that as A —2~ 
Io,2d diverges and the limiting value of J/t is two. Hence a critical value of Jc/t = 2 
is obtained in the 2d isotropic case, which is in agreement with the results obtained 
elsewhere [IHl [121 [S] • 



3.1. Three Dimensions - Isotropic Hopping 

The elliptic integrals of ( ITT]) become impossible to solve analytically in higher 
dimensions. However, to make progress with the three-dimensional problem it is 
necessary to understand as much about the properties of these integrals as possible. 
lo^nd is an odd function that is only defined for |y4| > and it is negative for A < —n. 
It is an important fact that the behaviour of /o,„d is qualitatively different for n < 3 
than for n > 3. This is because the limit as A ^ —n^ is divergent for n <3 but finite 
for n > 3, a result that can be obtained, e.g., by writing I^^nd in n-dimensional polar 
coordinates. 

To get a simple expression of J as a function of the energy A in three dimensions, 
analogous to f|T9|) . one can derive a set of equations analogous to f|T7j) . viz. 

=l{}t~ ^) 

IxX,3d = ^ AIj.^3^ — 2Ircy^3d 

Here we analyze the infinite 3d lattice with completely isotropic hopping frequencies 
and spin exchanges: J = J\\ = J± and t = t\\ = t± with = Ly = = oo. Starting 
from (|T2l) . and accounting for the symmetries in the elliptic integrals (J^. = ly = Iz, Ixy = 
lyz = Ixzi I XX = lyy = Izz ^ uotc that bclow, for simplicity, we refrain from using the 
3c? label for these integrals), one can eliminate C^^Cy and from the equation set to 
obtain the equation: 

(-2/,.,. J + 2h^JUh - QUilJ - AJhy + AJhyUh - 1 + Uh)Co = (21) 
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First we assume that Cq ^ and solve for U : 

2I,,JIo-6Ip + UI,yIo + Io ^ ' 

For finite A and J, the numerator of this expression is bounded - the / integrals are 
convergent for \A\ > 3 and finite. Thus taking the limit that [/ ^ oo is equivalent to 
setting the denominator to zero: 

2I,,JIo - 64V + 4JI,yIo + /o = (23) 

and then solving for J gives: 

J = -- ^2 (24) 

Finally, fl20|) can be used to give the simplest form for J, similar to f|T9|) . which again 
only depends on the elliptic integral Iq: 

J - ^^'^'^ (25) 
-1 + 4AIot ^^^^ 



The right-hand side of equation (1251) is a monotonically decreasing function as A 
approaches —3 from below, and therefore the minimal value of J is found from this 
limit. One can evaluate the right-hand side of this equation numerically to a very high 
accuracy. For example, one can complete the lattice sums of equation (fTTl) with larger 
and larger lattices in the limit of E approaching the bottom of the noninteracting band, 
that is -12t. However, a more accurate determination of this quantity can be obtained 
from the integral form of Iq, which, using the properties of elliptic integrals to simplify 
the multi-dimensional integral to an integration over a single variable, is given by 

27T (27r)2t(A + cosg) ^ ^ ^ 

where K{x) is the complete elliptic integral of the first kind. Again, examining the 
limit of A approaching -3 from below, and carefully accounting for the properties of the 
integrand near q = 0, one finds, in agreement with the lattice sum method mentioned 
above, that Io{A -3) = (-0.12633 ± 0.00009)/t Therefore, for this system we find 
JcA = 5.877 ±0.008. 

As mentioned, in the above derivation we have assumed for this solution that Cq ^ 
in writing down (!22l) . Putting Co = into the initial equation set, (fT2l) gives us either 
the trivial solution Co = = Cj, = = or a new equation for J. First we set 
Co = in (IT^ . and eliminating C^, Cy, from the equations yields 

J = ]^ (27) 

-A + AA'^Iot - 36hxt 

Numerical studies also indicate that this J has a similar qualitative behaviour as fl25|) . 
so if the limit as A ^ — 3~ of fl27|) is less than 5.877 then the initial calculation of 
Jc = (5.877 ± 0.008)t will be supplanted by the lower value. However, numerical 
calculations give a limit of J in ( !27|) of about 10.7. We note that this solution has 
a different symmetry than the first, a result that follows immediately from Co = 0; 
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so, this latter bound state is not an on-site or extended s-wave solution (see [H] for a 
discussion of analogous results in 2d). 

Finally, the critical spin exchange in the isotropic 3d case is determined to be 
Jc = (5.877 ± 0.008)t 7^ 2t. In part, we have provided this detailed analysis here 
because our results contradict a statement in [12] that the value of Jc/t = 2 (in the 
isotropic case) is independent of the lattice size and dimensionality. 

3.2. Two planes with no z hopping 

We now consider the anisotropic case of no z hopping in the lattice, for the reasons 
already reviewed in the introduction to this paper. Physically, one can see that to study 
an electron pair with no z hopping then one need only study the two plane problem since 
if the electrons are to see each other at all, they must be either on the same plane, or on 
adjacent z planes where they may interact via the spin exchange interaction Jy or J±. 
We only consider a single electron pair, and they are restricted to be in two neighbouring 
planes, so the problem collapses entirely to the two plane problem. Further, due to this 
restricted hopping the bottom of the two-electron scattering continuum will correspond 
to E = -8t. 

We start with the two plane problem, = Ly = oo and Lz = 2, that is isotropic 
in X and y and has zero z-hopping frequency. Here, as in f|T0|) - f|T2|) set = J. (The 
fact that for a spatial extent of = 2 it is necessary to specify open or periodic 
boundary conditions does not affect our results - see below.) Also, for the hopping we 
take t\\ = t and = 0. A great deal of simplification may be made with the / integrals 
of ffTTj) because the denominator of the integrand is independent of (and z is the only 
coordinate exhibiting anisotropy), and we can relate these integrals (which we denote 
by /) to their counterparts in the two-dimensional fully isotropic case. For example: 




mz=0,l qx,qy 



(28) 




Similarly 



-^0 — h,2d 

la = Ia,2d where a = x,y 

= Ia,i3,2d where a,P = x,y 
h =0 

laz = where a = x,y 
Izz = Io,2d 



(29) 



With these simplifications, (fT2l) simplifies significantly: 



Co — UIo,2dCo — 2JIx,2dCx — '^JIy,2dCy 



(30) 
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Cx = UIx^2dCo — '2JIxx,2dCx — '2JIxy,2dCy 

Cy = UIy^2dC0 — 2JIxy,2dCx — 2JIyy^2dCy 

Cz = — 2J±lQ^2dCz 

What is interesting about this set of equations is that the first three equations are 
identical to the 2d isotropic case, so for Cq 7^ we find that (fT9|) holds here as well. 
We can satisfy the third equation simply by setting Cz = 0. Since the z spin exchange 
is completely irrelevant to the equation, this physically means that this solution puts 
both electrons on the same plane, and thus the problem reduces completely to the 2d 
isotropic case. 

However, in our two-plane problem we can also obtain an equation for (in this case) 
J_L by allowing Cx = Cy = Cq = 0, and not get a trivial zero solution to the eigenvector 
of the Hamiltonian by having 7^ 0. Hence from fl30l) we immediately obtain: 

J± = -Ij^ (31) 

In this case Jy has dropped completely out of the equations, and this case corresponds 
to putting the two electrons on different planes, so they can only interact through 
J||. Recalling that Jo,2d diverges this gives us a critical spin coupling constant of 
Jc = J\\,c = 0. That is, we may obtain bound state solutions for arbitrarily small 
values J± (for any J). 



4. Discussion 



We have considered a model Hamiltonian motivated in part by the cuprate 
superconductors, and have examined the formation of two-electron bound states in 
the dilute limit. Work on related Hamiltonians fl2[ [T0| [Ti] and on other model 
Hamiltonians motivated by the cuprates p3] also followed this approach in the hope 
of better understanding the complicated many-electron states of these highly correlated 
systems. 

This paper analyzed the formation of electron pair bound states in a 3d lattice by 
comparing the minimum spin-exchange interaction that is necessary to allow for the 
formation of bound states for isotropic hopping and exchange, and for the complete 
suppression of interplanar hopping. We found that the minimum attractive interaction 
required is infinitesimal for the case where t± is zero, and hence the critical spin coupling 
for exchange-coupled planes is 

Jc = (32) 

The question remains of trying to extrapolate with this very simple model to the 
physics of two-electron bound state formation of the cuprate superconductors (albeit 
in the dilute limit for each plane). A 3d t ~ J — U model that takes into account the 
correct anisotropics of the hopping and superexchange, and for which t± is very small 
could indeed yield a two-electron bound state solution for very small values of J± if the 
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limit of Jc is continuous as ^ 0. However, in general J±^c will depend on both t± 
and J||, and the present study has not examined the parameter space associated with 
this much more difficult (numerically, that is) problem. (Note from equation ( 13T1) that 
there is no dependence on Jy when t±_ = 0, and for isotropic hopping and superexchange 
there is only dimensionless ratio, J/t.) 

Our results could be of interest to the study of ultra-cold atoms - see reference 
[16] and references therein for the connection of such systems to the kind of model 
Hamiltonians studied in connection with the cuprate superconductors. That is, in such 
systems one can tune the interactions, and with optical lattices the geometries, for 
such model Hamiltonians, so perhaps the limit that was examined in this paper can be 
created in such system. 

The above result necessarily leads to the question of the superconducting properties 
for such a model Hamiltonian in the dilute limit. As shown in [9], and as found for the 
case of exchange-coupled chains [17] , weak (or in our case zero) interplanar electronic 
hopping leads to a large density of states at half filling (the density of states for a 2d 
square lattice diverges logarithmically at the middle of band), necessarily leading to 
large increases in Tc when compared to the situation found when the hopping frequency 
is isotropic. Therefore, using a weak-coupling BCS approach one indeed expects for this 
model to show a large increase in relative to the 3d lattice. What happens in the 
strong coupling limit remains a topic for conjecture. In contrast to this result, note that 
for exchange-coupled chains [l?j the electronic density of states is in fact divergent at 
the band edges, and thus not surprisingly one finds the greatest enhancement of Tc at 
very low fillings. 
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